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ABSTRACT 

Context. RCW 120 is a well-studied, nearby Galactic Hn region with ongoing star formation in its surroundings. Previous work has shown that 

it displays a bubble morphology at mid-infrared wavelengths, and has a massive layer of collected neutral material seen at sub-mm wavelengths. 

Given the well-defined photo-dissociation region (PDR) boundary and collected layer, it is an excellent laboratory to study the "collect and 

■^C ■ collapse" process of triggered star formation. Using Herschel Space Observatory data at 100, 160, 250, 350, and 500 //m, in combination with 

r ^ ' Spitzer and APEX-LABOC A data, we can for the first time map the entire spectral energy distribution of an H n region at high angular resolution. 

* - , Aims. We seek a better understanding of RCW 120 and its local environment by analysing its dust temperature distribution. Additionally, we wish 

^ C"| to understand how the dust emissivity index, /3, is related to the dust temperature. 

O ! Methods. We determine dust temperatures in selected regions of the RCW 120 field by fitting their spectral energy distribution (SED), derived 
I ■ using aperture photometry. Additionally, we fit the SED extracted from a grid of positions to create a temperature map. 

Results. We find a gradient in dust temperature, ranging from > 30 K in the interior of RCW 120, to ~ 20 K for the material collected in the 

,£j , PDR, to ~ 10 K toward local infrared dark clouds and cold filaments. There is an additional, hotter (~ 100 K) component to the dust emission 

CZ2 . that we do not investigate here. Our results suggest that RCW 120 is in the process of destroying the PDR delineating its bubble morphology. The 

leaked radiation from its interior may influence the creation of the next generation of stars. We find support for an anti-correlation between the 

fitted temperature and/3, in rough agreement with what has been found previously. The extended wavelength coverage of the Herschel data greatly 

I , increases the reliability of this result. 
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1 . Introduction densations (plus an additional condensation, #9). They found ad- 

j 1 1 -■ ditional YSOs in the field, including a chain of 1 1 evenly-spaced 

RCW 120 d Rodgerset al. | | 1960|) is a Galactic Hn region that YSQs inside the most massive conden sation and a very dense 

O . displays a ring morphology at mid-infrared and sub-mm wave- core harboring a (poss ibly class 0) YSO. This work highlights 

lengths, and is presumably a bubble viewedm projection. It has the impact RCW UQ ig having on star formation far from the 

recently been studied by Zavagno_etaL| 420071 hereafter ZA07) ionizing source An analysis of the YSOs in the RCW UQ fieM 

and DeharvengetalJ (|2009, hereafter DE09) in the context of using the Herschel data pre sented here is given in a companion 

triggered star formation. It is only 1.3 kpc from the Sun (see pap er dZavagno et al.ll2010h . 
ZA07, and references therein), and is thus one of the closest 

Galactic H n regions. The ring morphology show n by RCW 120 i s a common 
ZA07 analysed the 1 .2 mm emission of RCW 120 and found feature of Galactic H n regions; | Churchwell et al ] ([2006, 2007) 
a fragmented layer of neutral material adjacent to the photo- identified almost 600 such ob]ects in the Sg/tor-GLIMPSE 
dissociation region (PDR). They identified eight millimeter con- data (Benjamin et al. 200j.| Deharveng et al. | d2010|) have shown 
densations, five of which lie in the collected layer of material, that over 85% of lnfrared < IR ) bubbles enclose Hn re g lons - 
but found no massive young stellar objects (YSOs) within the Because of their morphology, it is easy to locate the PDRs of 
condensations. They did, however, locate numerous YSOs sur- such bubbles and one can more easily identify the swept-up 
rounding RCW 120, indicating that star formation is active in material that is necessary for the collect and collapse process 
the region. Using APEX-LABOCA observations at 870/im, and ( |Elmegreen& Lada|1977|) . These objects present an opportunity 
Spitzer-MlPS observations at 24/an, DE09 extended this work to assess the efficiency of triggered star formation throughout the 
and calculated column densities and masses for the sub-mm con- entlre Galax y- Data from the Herschel telescope dPilbratt| | 2010D 
allow us for the first time to map the dust temperature variations 

* Herschel is an ESA space observatory with science instruments over an entire Hn region at high resolution. We can thus better 

provided by European-led Principal Investigator consortia and with im- determine the effect H n regions have on the creation of the next 

portant participation from NASA. generation of stars. 
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2. Data 

RCW 120 was observed by the Herschel telescope on 9 O ctober 
2009 with the PACS (Po ditschl I20l0h an d SPIRE dGriffinl 
2010) as part of, respectively, the HOBYS (Motte 201$ and 
"Evolution of Interstellar Dust" ( Abergel 2010) guaranteed time 
key programs. Data were taken in five wavelength bands: 100 
and 160 /mi for PACS (at resolutions of 10" and 14" and a final 
map size of 30' x 30'), and 250, 350, and 500 /mi for SPIRE (at 
resolutions of 18", 25", and 36" and a final map size of 22' x 
22'). We reduced these data using HIPE version 2.0. The SPIRE 
images used here are level 2 products, produced by the SPIRE 
photometer pipeline. We reduced the PACS data using a script 
provided by M. Sauvage. Recent calibration changes for the 
PACS and SPIRE data were taken into account (Swi nvardet al.l 
2010); throughout we use calibration uncertainties of 10% and 
20% respectively for the PACS 100 //m and 160 /mi bands, and 
15% for all three SPIRE bands dSwinvardet al.ll2010h. 

We also utilize Spitze r-MIPSGAL data dCarev et al.ll2009t) 
at 24 /mi and 70 /zm (at resolutions of 6" and 18"), and APEX- 
LABOCA data at 870/mi from D09 (at a resolution of 19"). 
The MIPSGAL 70//m maps have strong striping that we mediate 
with median filtering of the scan-legs (see Gordon et al. 2007). 
In Fig. Q] we show a 18' x 22' three-color image with SPIRE 
250 /mi (red), PACS 100 /mi (green), and MIPSGAL 24 /mi data 
(blue). Regions of interest are superimposed on the data; the con- 
densation numbers are from ZA07 and DE09. 



3. Dust properties 

We assume the dust emission in RCW 120 can be modeled by 
an optically thin grey-body, and that the emissivity of the dust 
grains can be fitted with a power law (Hildebrand 1983b : 



S v ocQB v (r) a-o|— N, 



(1) 



where S v is the measured flux density at frequency v, Q. is the ob- 
serving beam solid angle, B V (T) is the Planck function for tem- 
perature T at frequency v, Ko{v/vo~f is the dust opacity, ft is the 
dust emissivity index, and N is the dust column density. Each 
line of sight has variations in the dust temperature and emissiv- 
ity; the derived temperatures and emissivities therefore represent 
average values, weighted by the strongest emitting components. 
There has been much discussion about the value of the dust 
emissivity index, ft, and whether its value is anti-correlated with 
the dust temperature, T. Generally, ft is thought to range from 
1 to 2, but may also vary with wavelength (Meny et 



inge from 
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Dupac et al. (2003) found that ft and T are inversely related. 
Desert et a l. (2008) found a similar inverse-relatio n to exclude 
a consta nt value of/? at the 99.9% confidence level. IShettv etal.1 
(l2009allbr) . however, suggested that this relationship arises due to 
the influence of noise in a least-squares fit of Eq.Q] and from the 
combination of multiple emission components along the line of 
sight. A correlation between ft and T may indicate a change in 
dust properties at high density (Stenni k et al.ll2003l) . 

3. 1 . Aperture photometry 

To determine the dust temperature structure of RCW 120, we 
performed aperture photometry measurements on selected areas 
in the RCW 120 field. We resampled all images to the resolution 
of the MIPSGAL 24 /mi data to avoid pixel edge effects. There 
are three different areas of interest in the field of RCW 120: the 
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Fig. 1. Three-color image of RCW 120 composed of 250/mi 
(red) and 100 /mi (green) Herschel data, and 24 /mi MISPGAL 
data (blue). The entire field is 18' x 22'. Regions used for aper- 
ture photometry are shown and labeled. The coldest dust appears 
red, and is contained mostly in filaments on the southern edge of 
the bubble and within the infrared dark clouds to the north. The 
hottest dust appears blue and is largely confined to the interior 
region of the bubble, coincident with the ionized gas. The inte- 
rior of the bubble has little emission at wavelengths > 250 //m. 



interior of the bubble, the PDR (including condensations), and 
local infrared dark clouds (IRDCs) that appear in emission at the 
PACS and SPIRE wavelengths. The apertures we used are shown 
in Fig. Q] For each aperture, we selected a nearby background 
aperture that best characterizes the local background emission. 

For all apertures, we fit a single temperature to the coldest 
emission component of the spectral energy distribution (SED) 
using a non-linear least squares routine in two trials: once with 
ft — 2, and once with ft allowed to vary. For all apertures, the 
24 /mi emission is significantly greater than what would be pre- 
dicted by a single, cold temperature. This hotter component, 
likely caused in part by transiently heated small grains, is not 
well-constrained by the available data and therefore we ignore it 
in the present work. We also exclude the 70 /mi data point (and 
sometimes the 100 //m data point) if it is inconsistent with the 
emission from a single, cold temperature component. We leave 
the column density as a free parameter. 

The results of our aperture photometry are shown in TableQ] 
and examples of the SED fits are shown in Fig. [2] Uncertainties 
listed in Table[T]are the formal 1 cr values, taking into account un- 
certainties in the photometry measurements and calibration un- 
certainties. The direction of the interior of RCW 120 has dust 
temperatures of > 30 K. This result is rather uncertain, how- 
ever, because of the lack of emission above the background 
level in the interior of RCW 120 at wavelengths > 250 /mi. 
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Fig. 2. Example SED fits for the "Condensation 1" (red), 
"Interior 1" (blue) and "Condensation 5" (green) apertures. 
Points included in the SED fits are shown filled; points excluded 
from the fits are shown open. The fit with f3 allowed to vary is 
shown with a solid line and the fit with fi - 2 is shown with a 
dashed line. Apparent is the difference in temperature between 
the three apertures, and the higher temperature dust detected at 
24 pm (and 70pm and 100pm for Condensation 5). 

Table 1. Dust properties derived from aperture photometry 



Name 


T(K) 

/3 = 2 


T(K) I p 
ySfree 



Interior 


30.4 ± 0.9 


28.9 ±2.0 


2.4 ±0.2 


PDR - West 


27.1+0.6 


20.5 ±1.0 


2.9 ±0.2 


PDR - East 


23.1+0.7 


22.9 ±2.1 


2.1 ±0.2 


Cond. 1 


20.1+0.3 


21.6 ±1.1 


1.7 ±0.2 


Cond. 2 


22.2 ± 0.4 


22.1 ± 1.2 


2.1 ±0.2 


Cond. 3 


21.5+0.3 


22.2 ±1.2 


2.0 ± 0.2 


Cond. 4 


21.0 + 0.4 


23.0 ±1.2 


1.8 ±0.2 


Cond. 5 (IRDC) 


12.6 + 0.5 


9.2 ±0.9 


3.2 ±0.4 


Cond. 6 (IRDC) 


13.1+0.5 


10.5 ±1.2 


2.8 ± 0.4 


Cond. 7 


21.8 + 0.4 


22.7 ±1.2 


1.9 ±0.2 


Cond. 8 


14.3 + 0.6 


9.9 ±1.0 


3.3 ±0.4 


Cond. 9 (IRDC) 


12.8 + 0.5 


10.2 ±1.1 


2.9 ±0.4 


Cond. 10 


23.5+0.7 


25.9 ±2.8 


1.9 ±0.2 


Southern Filaments 1 


14.6 + 0.7 


10.7 ±1.3 


3.0 ±0.4 


Southern Filaments 2 


13.5+0.6 


9.6 ±1.0 


3.2 ±0.4 



The lack of emission at long wavelengths is evidence for a flat- 
ten ed, two dimensional ring instead of a bubble, as suggested 
bv iBeaumont & Williams (2010). The condensations and PDR 
have characteristic temperatures of ~ 20 K and the IRDCs have 
temperatures of ~ 10K, on the low end of the temperature dis- 
tribution for IRDCs found by Peretto (2010). The temperatures 
change little in the two trials. Uncertainties in temperature are 
small (often < 5%), and errors in/? are generally 10%. 

The (3 and T values, tog ethe r with the empirica lly derived 
curves of Dup ac et al.l (120031) and lDesert et alj (120 08). are shown 
in Fig. [3] Figure [3] clearly shows two T-fj groups, one for the 
PDR (and the condensations along the PDR) at (7\/3) - (20, 2.0), 
and one for dark clouds and cold filaments at (T,/3) =* (10, 3.0). 
We find an anti-correlation between/? and the dust temperature. 
Fitting a regression line of the same functional form as that of 
Desert et al. (2008) reveals: 
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Fig. 3. The relationship between the fitted temperature and dust 
emissivity index /?. The points are the results of our aperture pho- 
tometry measurements. The solid curve shows the best-fit line, 
the do tted curve shows the relationship found by Dupac et al. 
( 200 3J), and the da shed curve shows the relationship found by 
Deser t et all d2008l) . 



The f unctional form of the relationship given in iDupac et al.l 
(l2003h does not fit our data. While these data do support the 
f3 — T relationship, this result should be interpreted with caution 
because the calibration for the Herschel instruments is not yet 
finalized. Also, the fits to many of the cold filaments are con- 
sistent with a j3 value of two. For example, the j3 — 2 fit for 
Condensation #5 is also acceptable as it falls roughly within the 
formal errors (see Fig.|2]i. 

3.2. Temperature map 

While aperture photometry is well-suited for investigating the 
average dust properties of small areas in the field of RCW 120, 
a higher resolution map of dust temperature is useful for deter- 
mining the small-scale variations in dust properties. To construct 
such a temperature map, we fit the SED extracted at each pixel 
with a single-temperature grey-body model (Eq.[TJ, again leav- 
ing the column density as a free parameter. Sch nee et al.l {2007) 
showed that when attempting to simultaneously constrain f3 and 
T, uncertainties in the derived parameters can be very large due 
to the effect of noise in the data. We find that this is indeed the 
case and therefore in the following we fix /? to a value of two, 
which is valid for most areas of the map (see TableQ}. 

Fitting the SED of each pixel requires matching the resolu- 
tion and pixel locations of all wavelengths to those of the lowest 
resolution image. The SPIRE 500 pm data point has little impact 
on the derived fit parameters because we have higher-resolution 
data at 870 pm; we exclude the 500 pm data point so we can re- 
bin to the higher resolution of the SPIRE 350pm data instead. 
We smooth all images with a two-dimensional Gaussian rep- 
resenting the SPIRE beam at 350pm, and rebin to the 350pm 
resolution using MontageQ. We subtract an average background 
value from the data at all wavelengths, estimated from a field 
nearly devoid of emission. For all but the coldest dust, the Spitzer 
70 pm data point is dominated by emission from the cold com- 
ponent; we use this data point in the SED fits here. The inclusion 
of these data, however, causes the fitted temperature to be over- 



/3=(8.4±1.9)XT 



-0.4±0.1 



(2) 



http://montage.ipac.caltech.edu/ 
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estimated in very cold regions (see fits in Fig. [2]). 

The temperature map is shown in Fig. [4] for the same area 
shown in Fig. Q] We find the same results found with the aper- 
ture photometry method: the interior of RCW 120 is ~ 30 K the 
PDR is ~ 20 K, and the dark clouds are ~ 10 K. The highest 
temperatures (up to ~ 40 K) are found toward the southern edge 
of the bubble, surrounding the location of the main ionizing star. 
Errors in the derived temperatures for individual pixels are gen- 
erally < 20%. The exact location of the ionizing source is devoid 
of hot dust, an effect that may be caused by the action of stellar 
winds. We note, however, that the ionizing star of RCW 120 ha s 
been found to have relatively weak winds (Mart ins et alj |2010). 

The action of the ionized material appears to be creating 
openings in the PDR in several locations, shown with green lines 
in Fig. |U These openings can also be seen in Fig. [T] There ap- 
pears to be dust that is hotter by ~ 5 K outside the PDR at 
these locations, heated by radiation leaking through the holes 
in the PDR. This higher temperature dust may lead to the col- 
lapse of future generations of stars as it interacts with local 
molecular material outside the PDR. For the opening due south, 
this effect may have aided the creation of the sources inside 
Condensation #8 seen in Fig.[TJ the directions of the other open- 
ings show no indication of increased star-formation activity. 

There are a number of cold filaments in the field of RCW 120 
that have dust temperatures of ~ 10 K. All such cold filaments 
appear bright at wavelengths > 250 /im and therefore appear red 
in Fig. Q] These filaments are associated with IRDCs seen in 
Spitzer 8.0//m data. The most obvious cold patches in the field 
are from the IRDCs to the north, but there are numerous, thin 
filaments throughout the field. A number of cold filaments are 
oriented radially away from RCW 120, the coldest of which are 
towards the south-west. We suggest that these southern filaments 
are shaped by radiation from RCW 120 leaving the ionized re- 
gion. We suggest that the boundary of RCW 120 does not entrap 
all the emission; where the radiation leaks out, the filaments are 
compressed into the radial segments we see. In the direction of 
large condensations, the radiation is confined to the bubble area. 

4. Conclusions 

We have analysed the dust properties of the nearby bubble H n 
region RCW 120 using aperture photometry and temperature 
maps of Herschel, Spitzer, and APEX-LABOCA data. We have 
found a gradient of dust temperature, from ~ 10 K for local in- 
frared dark clouds, to ~ 20 K for the PDR (including sub-mm) 
condensations, to > 30 K for the interior of the bubble. 

Our results show support for a power-law form of the anti- 
correlation between the dust temperatu re T and t he dust emissiv- 
ity index/3, in good agreement with lDesert et alj (2008). Because 
of the range of temperatures found in a small spatial area, pre- 
sumably all at the same distance only 1.3 kpc from the Sun, 
RCW 120 is an ideal location to test for the ft - T relationship. 
With the wavelength coverage of Herschel, we are able to simul- 
taneously constrain f3 and T for regions of cold dust. 

The temperature map of the RCW 120 field reveals numer- 
ous locations in the PDR where radiation appears to be leak- 
ing into the surrounding interstellar medium. This implies that 
RCW 120 is in the process of fragmenting its nearly complete 
PDR layer. As the radiation leaves the confined environment, 
the radiation pressure shapes local cold filaments and may aid in 
collapsing local condensations outside the ionization front. 
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Fig. 4. The temperature map derived from the SED fits of each 
pixel, over the same area as Fig. Q] When deriving these maps, 
we fixed fi to a value of 2. The interior edge of the PDR, found 
by tracing the Spitzer 8.0/im data, is shown as the green line. 
The location of the ionizing star is shown with the small green 
circle. The location of hotter dust outside the ionization front is 
shown with the three green lines. 
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